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We propose a new method in the spectral analysis of noisy time-series data for damped oscillators. 

From the Jacobi three terms recursive relation for the denominators of the Fade Approximations 

built on the well-known Z-transform of an infinite time-series, we build an Hilbert space operator, a 

J-Operator, where each bound state (inside the unit circle in the complex plane) is simply associated 

^^ , to one damped oscillator while the continuous spectrum of the J-Operator, which lies on the unit 

("—^ ■ circle itself, is shown to represent the noise. 

f^ ' Signal and noise are thus clearly separated in the complex plane. 

r\l , For a finite time series of length 2N, the J-operator is replaced by a finite order J-Matrix Jjv, 

having N eigenvalues which are time reversal covariant. 

Different classes of input noise, such as blank (white and uniform), Gaussian and pink, are dis- 
cussed in detail, the J-Matrix formalism allowing us to efficiently calculate hundreds of poles of the 
Z-transform. Evidence of a universal behaviour in the final statistical distribution of the associated 
poles and zeros of the Z-transform is shown. In particular the poles and zeros tend, when the length 
of the time series goes to infinity, to a uniform angular distribution on the unit circle. Therefore at 
finite order, the roots of unity in the complex plane appear to be noise attractors. 

We show that the Z-transform presents the exceptional feature of allowing lossless undersampling 
r^ ' and how to make use of this property. 

C^ , A few basic examples are given to suggest the power of the proposed method. 
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I. INTRODUCTION 

Spectral analysis of highly noisy time-series data impacts such diverse fields as gravitational wave detection, Nuclear 
Magnetic Resonance Spectroscopy as applied to nuclear waste, brain/breast cancer detection, oil detection and other 
similar areas of application. 

Experimental time-series are always affected by the presence of noise. As long as the signal to noise ratio is not too 
poor, several filters are available to denoise the data within the framework of Fourier analysis and its variants. All 
such techniques, on the other hand, have drawbacks. 

Fourier analysis, being linear, treats noise on the same foot as the signal, and therefore is per se unable to distinguish 
the two. In the case of stationary uncorrelated noise, Wciner-Khintchine theorem can be used to denoise the data [l|], 
but it is unable to distinguish signal from correlated or non stationary noise. 

Wavelets denoising methods require some knowledge of the signal to be found in order to avoid smoothing out the 
signal itself \2\ . 

Prony's method is tailored to analyze signals composed of damped oscillators, but it is known to be unstable in 
the presence of noise. This is because -at least in its original formulation- it assumes an "all-pole" system; in other 
words, it is equivalent to the construction of a Fade approximant (rational approximation) with constant numerator. 
Having no zeros it is unable to model noise, as we shall see in section Hi B 31 Modifications have been proposed to 
stabilize it Q, but, apart from being cumbersome, no effort has been made to classify the poles so as to distinguish 
the noise ones from the signal ones. 

The matched filtering technique, which calculates cross-correlations between noisy detector outputs and reference 
waveforms from a "library" , is only useful when the waveforms are well predicted theoretically [4] . 

Still, all the above techniques fail when the (average power) signal to noise ratio approaches 1. It is exactly this very 
high noise case we deal with here, presenting a denoising method based on the analytic properties of the Z-transform 
(generating function) of the data Q . 

The present paper is thus organized: 

From a given infinite time series, we build -by way of Fade Approximations of its Z-transform- a tridiagonal Hilbert 
space operator, a J-Operator, which, for a finite time series of length 2N, is replaced by a finite order J-Matrix 3n, 
having N time reversal covariant eigenvalues. Such eigenvalues correspond to the poles of the Z-transform itself. 

We then introduce a specific class of signals made of an arbitrary number of damped oscillators. The poles of the 
Z-transform of such signals are known to be inside the unit circle [5| . 

It is on the other hand known that a Taylor series with random coefficients has the unit circle as natural boundary 

t3- . . . . 

For a time series of noisy damped oscillators, the spectrum of our J-Operator will therefore comprise a discrete 

spectrum where each bound state (inside the unit circle) is simply associated to one damped oscillator; and a continuous 
spectrum which represents the noise, and lies on the unit circle in the complex plane. 

It is known that for a finite time series the natural boundary due to the noise is replaced by doublets of poles and 
zeros [Ij, [ij, [l^, [lg| surrounding the vicinity of the unit circle. We numerically study the statistical distribution of 
such poles and zeros for purely noisy time-series, considering different input classes, such as blank (white and uniform), 
Gaussian and even coloured noise. We find compelling evidence of universal behaviour. In particular the poles and 
zeros tend, when the length of the time series goes to infinity, to the roots of unity in the complex plane which appear 
to be noise attractors. 

As the Fourier Transform can be seen as the specialization of the Z-transform to the unit circle, the above result, 
combined with the knowledge that the poles associated with damped oscillators are inside the unit cicle, explains the 
poor behaviour of the Fourier Transform and the like while analyzing damped oscillators spectral data corrupted by 
heavy noise. 

We then prove that the Z-transform presents the exceptional feature of allowing lossless undersampling. This 
allows us to use "interlaced sampling" to obtain multiple time sequences from a single one and thus improve detection 
through coincidence. 

Finally, we give some simple examples to suggest the power of this approach. As is usually done in Nuclear Magnetic 
Resonance Spectroscopy, we consider one or more damped oscillator signals repeatedly observed. By covering the 
complex plane with a lattice and counting the "non-zero paired" poles in each cell (that is poles which have no 
associated zero in their vicinity), an approximate evaluation of both the detection probability and the false alarm 
level for a single observation can be obtained. "Interlaced sampling" can, on the other hand, significantly improve 
performance even in the case of a single time sequence. 



II. ANALYTIC PROPERTIES OF THE Z-TRANSFORM 

Given a time-series sq, si, S2, Sn, , it is a standard practice ^5,] to associate to it its Z-transform defined as: 



Z(z) = ^s„z-" 



(1) 



Clearly, due to the boundedness of the time-series coefficients, this function is at least analytic outside the unit 
disk. 

A. Construction of the tridiagonal Jacobi Matrix. 

We shall now introduce the Hilbert space operator whose resolvant matrix element is precisely the Z-transform. 
We are going to deal mainly with the [^^^] (z) Pade Approximant to the Z — transform of a given time series 
and write 



n ~ 1 



(^) 



Nn-ljz) 
Dn{z) 



(2) 



Where Dn{z) is a monic polynomial of degree n and A^„„i(z) a polynomial of degree n — 1 in z. 
It is known [7] that Dn{z) satisfy the Jacobi three terms recursive relation 



where 



Dk+i{z)^[z~Ak]Dk{z)^RkDk-iiz) D^i{z) = l Doiz) ^ I 



Ak = -{a2k + a2k+i) Rk = a2k-ia2k fc > 1 oq = 



(3) 



(4) 



The ak being the coefficients of the Stieltjes continued fraction expansion of the Z — transform . 
This relation can be written as an eigenvalue problem by introducing the tridiagonal {N + 1) x [N + 1) Jacobi 
Matrix 



Jn — 



Then, when Dn^i{z) = 0, the previous recursive relation can be written as 
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(5) 



where V is the column vector 



Jj^V = zV 



V'^ = [Da{z),Di{z),...DN{z),] 



(6) 



(7) 



We thus see that the zeros of Dn+i{z) are the eigenvalues (zq, zi, ..., zn) of the tridiagonal Jacobi Matrix of order 
N+1. Therefore, the Z — transform appears to be a matrix element of the resolvant of the J — Matrix. 
We also introduce the norm of the eigenvector: 



p=N 
p=0 



(8) 



B. The case of infinite time-series 

1. Damped oscillating signals Z-transform, 

One of the most important application of the Z-transform concerns damped oscillating signals, such as those found 
in nuclear magnetic resonance data [8] and in burst gravitational waves [9|, |lO| . 

Let us start by considering the signal in the absence of noise: for a finite ensemble of damped oscillators, the 
discretized data will read 



Sk 



^^Ape''^"^^ fc = 0,l,2,...,7V-l ujp = 2t: fp + iup (9) 



where Ap, fp, and ap are the amplitude, the frequency and the damping factor of the p oscillator. T is the recording 
time and N the number of data. 

We now consider their Z-transform ^ and allow both the recording time T and the number of data N to go to 
infinity, while keeping T/N constant, so that ^ becomes the (formal) Taylor series : 



+ 00 



Z(z) = ^SfcZ-'= (10) 



fe=0 

Combining the previous two equations, we get for the Z-transform the expression: 



+ 00 +00 +OC . 

z{z) = y: ^^--' = E -"' E ^ve'-'^^ = E ApY^iz^'e-^^' = Y. 7-44:;;^ (^i) 

fc=o fc=o p p k=o p ^ e "T^ 

From this we immediately see that, in this noiseless case, the Z-transform is a rational fraction in z, the poles of 
which are 

Zp = e''^"^ (12) 

Notice that all poles have to lie strictly inside the unit circle because Imc^p > 0. The residues are 

Pp = ZpAp (Ij) 

Therefore, after identifying the poles Zp and residues pp of the Z-transform, we get the quantities of interest, which 
are the frequencies 

the damping factors 



and the amplitudes 



N 



A, = ^ (16) 

Zv 



of the damped oscillators of the signal. 

This straightforward calculation shows that for a sum of oscillating damped signals, the Z-transform associated to 
its time-series is a sum of poles in the complex plane. The position of each pole is simply linked to the damping factor 
and the frequency of each of the oscillators. Also, it is important to note that all these poles lie strictly inside the 
unit disk. 



2. Noise Z-transform 

Due to the linearity of the Z-transform, the Z-transform of a noisy time-series will be the sum of the Z-transform 
of the signal plus the Z-transform of the noise. 

It is therefore of the uttermost importance to get a deep insight inside the analytic structure of a purely noisy 
time-series. 

The following theorem is the key to the answer to this question. 

Steinhaus theorem: A Taylor series with random coefficients has, with probability one (that is: except for a set 
of measure zero), the unit circle as natural boundary [11| . 

3. Tfie case of noisy data 

Combining the previous statements, we get the following simple but hitherto unrecognized fundamental result. 

The Z-transform of a finite number of damped oscillations is, in the presence of noise, an analytic function outside 
the unit disk, plus a rational function having a finite number of poles strictly inside of the unit disk, each of them 
representing one damped oscillation. 

We can state the following fundamental theorem. 

Theorem (Bessis-Perotti). 

To any noisy time series can be associated a tridiagonal Hilbert space operator, its J — Operator, extension when 
A'^ — > oo of the previously defined J^r. With probability one (that is: except for a set of input data series of measure 
zero), the spectrum of this operator is made of two parts: 

An essential spectrum with support the unit circle. This spectrum is associated with the noise (uncorrelated part 
of the signal). The corresponding eigenfunctions have infinite norm. At a finite order, when the infinite J — Matrix 
is truncated, it decomposes into the Froissart poles of the Froissart doublets which will be described in the following 
section. 

A discrete spectrum, made of a finite number of poles inside the unit circle: each pole represents a component of 
the signal made of a finite number of damped oscillators. The corresponding eigenfunctions have finite norm. 

C. The case of finite time-series 

Up to now we have been dealing with time-series of infinite length, we come now to the more realistic case where 
the time-series has a finite length N. 

The problem of constructing the Z-transform Z(z) of such finite time-series is a standard problem in Mathematics 
known as the Fade Approximation construction ^] . 

1. The special case of non noisy finite time-series 

If P is the number of damped oscillators, as soon the length of the time-series is A^ > 2P + 1, the diagonal Fade 
approximant P/P gives an elegant solution to the above problem. 

2. The general case 

In this case, M. Froissart [12| has shown that the natural boundary generated by the noise is approximated by 
doublets of poles and zeros (Froissart doublets - see e. g. [l^ [ij, ^^ [lg| ) surrounding the vicinity of the unit circle. 
The mean distance of the members of each doublet being of the order of the noise magnitude and we conjecture that 
it decreases exponentially with the length of the signal. 

III. UNIVERSALITY PROPERTIES OF THE NOISE PART OF THE SIGNAL 

Here, we are mainly interested in studying the statistical distribution (both radial and angular) of the Froissart 
doublets as a function of the input noise statistical properties and the length N of the time-series. 



The most unexpected result is the universahty of the distribution: the radial distribution in the purely noisy case 
is universally Lorcntzian, the phase distribution is uniform (approaching the [Af/2] roots of unity when N goes to 
infinity) . 

Due to this property of universality, the only parameter left is the width of the Lorcntzian, which depends on the 
length N of the time-series. A conjecture, based on heuristic arguments, about the functional dependence of this 
width on N will be provided in the sequel. 

A. Experimental evidence for universality 

We calculated the radial distribution of the diagonal and sub diagonal Pade approximants for time-series of up to 
600 data, which allow the construction of Pade approximants up to 300/300. 
To test universality, we have been studying 7 different kinds of noise: 

1) the complex input data have uniform distribution both in modulus and phase. 

2) the complex input data have uniform distribution inside a square domain. 

3) the complex input data have uniform distribution inside a circle domain. 

4) standard pink noise. 

5) the complex input data have Gaussian distributions. 

6) the complex input data are correlated following an autoregressive moving average (ARMA) model. 

7) the input data are real and have a Gaussian distribution. 

Typical results for the single run radial distribution of the Pade approximant poles in some of the above seven cases 
are shown in Figure ([T]): in all cases the behaviour is Lorentzian and centered on 1. 

Figure ([2|) presents instead the phase distributions of the same Froissart poles: in all cases the poles are uniformly 
distributed between and 27r, the deviation being Gaussian. 

The graphs shown refer to time-series of 200 data: the case for which we have run the highest number of tests, 
varying several times the seeds of the random number generators used for each of the 7 kinds of noise. 

Analogous results apply to the zeros of the Pade approximants, with the difference that the radial distribution is 
centered on a point larger than 1 which approaches unity when N goes to infinity (at present, we have no explanation 
for this slightly unsymmetrical behaviour between the poles and zeros statistical distributions). 

We have recently gained evidence that Universality extends to the case of noise statistical distributions having 
infinite second moment and, even, infinite mean. We shall report the details in a forthcoming paper. 

For the very special case of Gaussian uniform noise, a mathematical proof of our results can be found in Ref. [l7|. 

B. The width of the universal radial Lorentz distribution 

We are left with the problem of finding an expression for the width of the universal radial distribution of both the 
poles and zeros and for the distance of the center of the radial distribution of the zeros as a function of the time-series 
length N. 

On heuristic arguments, we propose for all three cases the following formula 

VK.oM^ (17, 

Figure ([3]) shows a fit to the above formula of said three quantities. The agreement appears quite good and seems 
to improve when increasing TV. 

IV. FIRST CONCLUSIONS 

Spectral analysis has been dominated by the Fourier transform (and its variants), which is the ideal tool when noise 
is not too important. 

The major drawback of the Fourier transform is its linearity which prevents it to distinguish noise from signal 
because the Fourier Transform treats the noise on the same foot as the signal, and therefore leaves noise intact. 

Furthermore it is clear that the discrete Fourier Transform is nothing but the specialization of the Z-transform to 
specific points of the complex plane which are the roots of unity. This is dramatic because, as we have shown in our 
presentation, the roots of unity are the noise attractors of the complex plane. 

This explains the poor results of data spectral analysis in the presence of heavy noise produced by the FFT. 



We therefore suggest to substitute the Z-transform for the discrete Fourier Transform when analyzing highly cor- 
rupted time-series data. 

V. BEYOND THE FOURIER ANALYSIS: METHODS FOR EXTRACTING SPECTRA FROM NOISY 

DATA. 

In order to make a clear distinction among the poles of the Z-transform of a noisy time series of length 2N + 2, 
and identify the eigenvalues corresponding to the discrete spectrum, we can make use of several properties of the 
Z-transform. Each method of those we are going to present addressing a different property, a comparison of the 
results of all methods applicable to any specific time-series is likely to be the optimal choce. 

For some of the methods delineated below, we need to find both poles and zeros of the Z-transform. For the former 
ones, we compute the {N -|- 1) x (A^-|- 1) Jacobi Matrix associated to the denominator of the Z-transform we described 
in section Hi Al and then diagonalize it. The latter ones can be calculated diagonalizing the N x N Jacobi Matrix 
for the numerator of the Z- transform Fade Approximant, which is simply the denominator Jacobi Matrix with the 
first column and first row removed. The procedure for computing both zeros and poles of the Z-transform therefore 
consists in the diagonalization of an already tridiagonal matrix; this makes it computationally extremely efficient, and 
allows us to easily calculate hundeds of poles and zeros. 

A. Cleaning of Froissart doublets 

A first denoising method consists in removing from the complex plane the poles which can be identified as part of 
Froissart doublets, due to their proximity to zeros of the Z-transform. 

This can be done by ordering the poles and zeros in couples in order of increasing distance between the two, and 
keeping only those whose distance is higher than a given value e (usually e > 0.2). 

B. Variational principle 

We can on the other hand, making use of the results of section III Al search the eigenvalues of Jacobi Matrix 
associated to the denominator of the Z-transform having modulus strictly smaller than one and for which the norm 
of the eigenvector have local minima. Such eigenvalues correspond to the frequencies of the signal. 

C. Stationariety 

It is also possible to compare Fade approximants of different order and look for stable "non-zero paired" poles 
inside the unit circle: these will be the signal poles, while the non stationary poles will be linked to the noise. This 
method is particularly useful in dealing with poles in the tails of the noise pole distribution. 

D. Undersampling 

Another method makes use of an exceptional feature, never reported up to now, of Fade Approximations with 
respect to undersampling. 

Emmanuel Candes |18l | has extended the Nyquist-Shannon sampling theorem [19] and shown that it was possible 
under some circumstances to use bandwidth smaller than the ones requested by the Nyquist theorem (sufficient 
condition only but not necessary) without losing information. 

It is therefore interesting to look for classes of signals for which the corresponding Z-transforms allow undersam- 
pling without any loss of information. We are going to show that the class of damped oscillators allows lossless 
undersampling. 

Suppose that in (jlOp . we undersample the time series by taking only the Sk with k — mn where m is given fixed, 
while n runs from zero to infinity to start with. We get the undersampled Z,„ — transform 

+ OC 

Z™(z) = ^s™„0-" (18) 

n=0 



hoc +00 +00 



^™(.) ^ E^™---" - E -"" E -^.^'^-^" - E -^^ Y.i^-^e-^^r = E , ,_t„.,x (19) 

n=0 n=0 p p n=0 p -*- "^ '^ 

We see that the Pade Approximation P/P again gives the exact solution to the undersampled time series. So 
provided we pick N — 2P + 1 equispatiated terms in the time series, we are able to rebuild it exactly. 

We point out here that we assume to be dealing with data series such that the periods of the signal oscillations to 
be reconstructed cover several data points. This means that, as long as m is less than the number of data points per 
period of each of the signal oscillators, the reconstructed signal frequencies will be given by the mth roots of the poles 
Zp having the lowest phases. 

The above result suggests another technique -which we call "interlaced sampling"- to improve sensitivity when only 
a single data time-sequence is available. Taking advantage of the undersampling properties of the Z — transform, 
we divide the data in m undersampled sequences, the first one comprising the points 1, to -1- 1, 2m + 1, 3to -|- 1, ..., the 
second one comprising the points 2,m + 2, 2m + 2, 3m -I- 2, ... and and so on. To first oder, the ratio of the width of 
the pole distribution around the unit circle to the distance of the signal pole from it does not change, but instead of a 
single signal pole, we do now get a cluster of m poles, which are much easier to detect. Preliminary tests have given 
encouraging results. 

VI. APPLICATION OF THE METHOD TO EXPLICIT EXAMPLES 

We shall now proceed to describe some simple examples. The actual procedure we use to reconstruct the signal is 
very basic, but sufficient for our present aim of showing the potentiality of our proposed system. Plain Fortran was 
used to write the code implementing the procedure described in the following. 

A. A first example 

As a first example, we consider a damped signal at zero frequency. The chosen damping factor reduces the signal 
by a factor 10^^ over 300 signal points, resulting in a pole at (0.95238,0.0). Noise is of the type 1) above (the complex 
input data have uniform distribution both in modulus and phase). 

A signal series of 300 points is relatively short. Since the separation of signal and noise is only univocal in the limit 
of an infinite series, and -moreover- a single pole off the unit circle could be due to the particular noise sequence used, 
we repeat the sequence several times, as is for example usually done in Nuclear Magnetic Resonance Spectroscopy. In 
the present case, we repeat the sequence 40 times, each time changing the seeds for our noise generator, and search 
for a stable pole inside the unit circle with no zero close to it. 

For each of the 40 data sequences we clean most of the Froissart doublets according the procedure delineated above: 
we calculate the distance of each pole from all the zeros, arrange the pole-zeros couples in order of increasing distance 
and keep only the uncoupled pole and the poles whose distance from the coupled zeros is higher than a set value. 
We find that cut distances between 0.6 and 0.8 -which leaves us with two or three "non-zero paired" poles per run- 
give good results for the range of signal/noise ratios we explored; for higher signal/noise ratios good results can be 
obtained with much smaller cut distances. 

Directly plot of these poles (and their coupled zeros) for all 40 runs, can be misleading to the eye. To better see 
areas of concentration of poles, we divide the complex plane in a number of square boxes and plot in three dimensions 
the number of poles minus the number of coupled zeros in each box. For the tests we present, the boxes were of side 
0.1. 

Figures ^ and ^ show the case of an average power signal/noise ratio (SNR) of approximately —lOdb (matched 
filtering SNR equal to 5.5), using a cut distance of -respectively- 0.75 and 0.2: in both cases the highest peak is at 
the position of the signal pole. Note that this can be considered a worst case scenario: the signal pole lies close to a 
"four corners point". The reconstructed poles therefore are distributed over four boxes: for a cut distance of 0.75 the 
four peaks sum to 17, while the height of the highest "fake" peak is only 2. For a cut distance of 0.2 the four peaks 
sum to 22, but the height of the highest "fake" peak is now 6. This suggests that reducing the cut distance this much 
is not convenient: going from 0.75 to 0.2 improves the detection probability from about 40% to about 50%, but at 
the same time forces us to increase the false alarm threshold by a factor 3 (from 2 to 6 out of 40). For the examples 
that follow we shall therefore show only the results obtained using 0.75 as our cut distance. 

Reducing the average signal/noise ratio to approximately —15db (matched filtering SNR equal to 3.0), the peak is 
still the highest, but "fake" peaks far from it become more visible, as can be seen in Figure ©. This is mostly due 



to a higher noise induced spread of the reconstructed poles: if we sum the four signal peaks we get 15, which is not 
much less than the 17 of the previous case (detection probability is still around 40%), especially considering that the 
height of the highest "fake" peak is still 2 which means we can keep the same false alarm threshold. 

Figure ([7]) finally shows the case of an average signal/noise ratio of approximately —18db (matched filtering SNR 
equal to 2.2): the peak at the signal pole is still the highest, but peaks far from it (the height of the highest still being 
2) are now substantial compared to it and to the sum of the four signal peaks which is now only 9. The false alarm 
theshold does not need to be changed, but detection probability is now down to 25%. 

B. A more complicated signal 

Adding to the —15db case an oscillating term, resulting in a couple of complex conjugate poles in (0.67344, 0.67344) 
and (0.67344, —0.67344) respectively, and having a average signal/noise ratio of approximately —21db (matched filter- 
ing SNR equal to 1.5), gives us Figure ([8]). The oscillating signal, even if weaker, appears to give peaks much sharper 
than the peak given by the non-oscillating signal, but only because its poles are far from box borders: if we sum the 
four peaks corresponding to the non-oscillating signal we get 19, the same as the highest oscillating signal peak (the 
other one is only 15, the asymmetry being due to the noise being complex). Note that the height of the highest "fake" 
peak is now 4, which suggests the need for a higher false alarm threshold. 

VII. GENERAL CONCLUSION. 

There is still plenty of room for improvement: of the methods outlined in section |V] we have used, in the examples 
given above, only the cleaning of Froissart doublets and the knowledge that the signal poles must be inside the unit 
circle. 

We are moreover still not using all the information available. For example we did not use the residues of the poles, 
from which the signal amplitude can be obtained, and we haven't yet performed a complete optimization over the 
various parameters used in the analysis. In particular there is the need for an extensive mapping of the spread of 
the reconstructed poles as a function of noise level and Fade approximant order. Our last example also suggests that 
there might be a connection between number of signals and fake alarm threshold, which needs to be investigated too. 

This said, we can nonetheless conclude the following. 

In our method, the Z-transform appears as an extension of the discrete Fourier transform to complex values of the 
frequency: in the complex plane the discrete Fourier transform is the restriction of the Z transform to values of z 
located at the roots of unity. However, the noise is not uniformly distributed in the complex z-plane and, unfortunately, 
the roots of unity are noise attractors as shown at the beginning of the article. Therefore, the Fourier transform is 
not a good choice when the signal is embedded in high noise as it is the case in many circumstances. The analytic 
treatment of the noise we propose, distinguishes, in a drastic way, the signal from the noise by their totally different 
analytic properties. The fact that the noise presents an analytic characterization that is universal, independent of its 
statistical properties, (the most novel result of our paper) has as consequence a much greater independence of the 
signal identification on the noise level. In principle, in the limit of very long signals, the signal extraction becomes 
independent of the signal to noise ratio. 

The examples we have given, moreover show that the Z — transform analysis promises not only in theory to be 
an extremely powerful tool especially for the frequently encountered case of damped signals in a highly noisy 
environment, as it gives us a sure signature for the signal by radially separating it in the complex plane from the 
noise. 
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FIG. 1: Radial cumulative distribution of the Froissart poles of the [100/100] Fade approxiniant for different kinds of noise (see 
text for a full description). Full line, black: type 1; dash, red: type 2; dot, blue: type 3; dash-dot, green: type 4; dash-dot-dot, 
magenta and dark yellow: two different cases of type 5. The thick full line is a fitting by a Lorentzian cumulative distribution 
((I/tt) arctan[(p — 1)/W] -f 1/2, where W is the width of he distribution). 




pole number 



FIG. 2: Phase distribution of the Froissart poles of the [100/100] Fade approximant for different kinds of noise (see text for a 
full description). Full line, black: type 1; dash, red: type 2; dot, blue: type 3; dash-dot, green: type 4; dash-dot-dot, magenta 
and dark yellow: two different cases of type 5. 
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FIG. 3: The width of the universal radial distribution of both the poles (squares, black) and zeros (up triangles, dark yellow), 
and the distance from unity of the center of the radial distribution of the zeros(down triangles, blue) as a function of the time- 
series length N. Also shown is their fitting to ll7l The global multiplicative factor a is, respectively, 0.404 ± 0.005, 0.52 ± 0.01, 
and 0.62 ±0.01. 
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FIG. 4: Signal/noise ratio of approximately — 10d6. Distribution of the [149/150] Fade approximant poles distant from their 
coupled zeros more than 0.75. The vertical axis represents the number of poles in each box of side 0.1. The highest peak 
corresponds to the signal pole at z — (0.95238, 0.0). 




13 



1.0 x^-^"^ 



1.0 
0.5 A 

,9'^ 



FIG. 5: Same as Figure |4] above but the minimum pole- zero distance has been reduced to 0.2. 
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FIG. 6: Signal/noise ratio of approximately — 15db. Distribution of the [149/150] Fade approximant poles distant from their 
coupled zeros more than 0.75. The vertical axis represents the number of poles in each box of side 0.1. The highest peak 
corresponds to the signal pole ai z — (0.95238, 0.0). 
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FIG. 7: Signal/ noise ratio of approximately — 18d6. Distribution of the [149/150] Pade approximant poles distant from their 
coupled zeros more than 0.75. The vertical axis represents the number of poles in each box of side 0.1. The highest peak 
corresponds to the signal pole at z — (0.95238, 0.0). 
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FIG. 8: Signal/noise ratio of approximately — 15d6 for the constant term and of approximately —21db for the oscillating 
one. Distribution of the [149/150] Pade approximant poles distant from their coupled zeros more than 0.2. The level curves 
represent the number of poles in each box of side 0.1. The two highest peaks correspond to the oscillating signal poles at 
z = (0.67344, 0.67344) and z = (0.67344, -0.67344)). The non-oscillating signal pole at z = (0.95238, 0.0) appears to be lower 
and more spread out. 



